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Abstract: "Brane Induced Gravity" is regarded as a promising framework for addressing 
the cosmological constant problem, but it also suffers from a ghost instability for param- 
eter values that make it phenomenologically viable. We carry out a detailed analysis of 
codimension > 2 models employing gauge invariant variables in a flat background approx- 
imation. It is argued that using instead a curved background sourced by the brane would 
not resolve the ghost issue, unless a very specific condition is satisfied (if satisfiable at all). 
As for other properties of the model, from an explicit analysis of the 4-dimensional graviton 
propagator we extract a mass, a decay width and a momentum dependent modification of 
the gravitational coupling for the spin 2 mode. In the flat space approximation, the mass of 
the problematic spin ghost is instrumental in filtering out a brane cosmological constant. 
The mass replaces a background curvature that would have had the same function. The 
optical theorem is used to demonstrate the suppression of graviton leakage into the un- 
compactified bulk. Then, we derive the 4-dimensional effective action for gravity and show 
that general covariance is spontaneously broken by the bulk-brane setup. This provides a 
natural realization of the gravitational Higgs mechanism. We also show that the addition 
of extrinsic curvature dependent terms has no bearing on linearized brane gravity. 
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1. Introduction and Results 

"Brane induced gravity" (BIG)(l|, §] (for a review see Q) is regarded as providing one of 
the few promising frameworks for addressing the cosmological constant problem [Q] in the 
sense of explaining how the observed value of vacuum energy could be so small as compared 
to the types of values inferred from quantum field theory j|, || (other approaches include 
Q. H|, for reviews see |], 1C, |ll|). In its basic form, the model regards the 4-dimensional 



universe as a 3-brane in a non-compact (4 + n)-dimensional bulk spacetime. Gravitational 
dynamics arise from Einstein-Hilbert actions both in the bulk and on the brane. Bulk 
gravity absorbs the effect of brane cosmological constant, leaving behind effectively 4- 
dimensional gravity on the brane with a "filtered out" cosmological constant. Unfortunately 
consistent constructions are still lacking, for example, issues related to the presence of a 
tachyonic ghost mode have not been resolved satisfactorily. Hence more work is needed to 
obtain consistent working models of this type 1 . 

1.1 An overview of BIG, its promise and problems 

Let's briefly review the origin, workings and shortcomings of the model. 

Historical background and motivation: In ordinary Einstein-Hilbert gravity a 
cosmological constant A in a 3 + 1 dimensional spacetime gives rise to a de Sitter metric. 
However, if this specetime is a 3-brane embedded in a 4 + n dimensional bulk, then it 
is known that bulk gravity has classical solutions in which A (which is the same as the 
brane tension) curves directions transverse to the brane while the metric along the 3-brane 
remains flat, provided n > 1 [fH| , |i"5| , |il| . This eliminates A from the brane point of view 
at the level of classical solutions. But in this model, gravity is a (4 + ra)-dimensional force. 

To get effectively 4-dimensional gravity on the brane without compactifying the extra 
dimensions, Dvali, Gabadadze and Porrati [II], 0j proposed adding to the bulk gravity action, 
a 4-dimensional Einstein-Hilbert action for the brane induced metric. Such brane localized 
actions, along with extra corrections, naturally arise in field theory |17| as well as in string 



theory [18, |19|, 20, 21, |22|, 23] setups, hence the name brane induced gravity. However, for 
phenomenological reasons, the parameters have to be tuned to make the brane localized 
action dominant. The n = 1 case is the well studied DGP model |l], |3|. It is special in that 
it does not suffer from tachyon/ghost problems. But it is also not adequate for addressing 
the cosmological constant problem in the sense stated above. In || |B| it was argued that 
models with n > 2 should be able to filter out a brane cosmological constant, making 
gravity on the brane less sensitive to it, while at the same time giving rise to essentially 4 
dimensional gravitational interactions on the brane. These arguments were mainly based 
on the classical solutions of [14, 15 1 described above 2 . 



Beside the classical solutions, most of the explicit work on BIG with n > 2 has fo- 
cused, for reasons of technical simplicity, on the A = case, analyzing the propagation of 



1 In this paper we will not discuss variants of the basic setup like cascading gravity WM that address the 
ghost issue p3[. 

2 The n — 2 models also filter out A but they share features of both n — 1 and n > 2 models and are 
discussed in some detail in |24, M, UM. 



linearized metric fluctuations for a flat brane, in flat bulk spacetime. The purpose is to, 
from this, glean information on infrared modifications of gravity that would eventually re- 
sult in filtering out A. In particular, one finds that for appropriate values of parameters in 
the action, gravity can be made effectively 4-dimensional over a desired range of distances. 
It will also contain massive unstable gravitons g, [5], || (27J. Beyond this range (say, for 
distances of the order of the size of the Universe), the emergent unstable graviton, with a 
small but momentum dependent mass, produces infrared modifications of gravity. This is 
supposed to be relevant to filtering out A, as suggested by the classical solutions mentioned 
above. 

Problematic issues: The above picture looks promising, except for some important 
issues that arise in the analysis carried out so far, and that need to be better understood 
and resolved. We describe five issues below. 

1. The evidence in H, E| that BIG would filter out a brane cosmological constant is 



mostly based on the existence of the classical solutions of [14, 15] that are supposed 
to describe how A curves the space transverse to the brane. However a closer ex- 
amination shows that in the vicinity of the brane, the source structure implied by 
these solutions does not correspond to the brane source (which is a common prob- 
lem of brane solutions). This means that a detailed analysis of fluctuations (needed 
to obtain the effective 4-dimensional gravity on the brane) cannot be carried out 
around existing classical solutions 3 . To avoid this problem, in this paper we start 
with a flat background approximation and try to glean information about the effects 
of background curvature sourced by the brane. 



Second, as pointed out in [28], the analysis of BIG in a flat background approximation 
indicates the presence of a tachyonic ghost mode and, therefore, an inconsistency of 
the model. This is a major drawback and resolving the ghost problem is crucial for 
developing these models further. It is easy to see that a ghost generically appears in 
massive gravity theories unless it is given an infinite mass, as in Fierz-Pauli gravity 



p9| , 30, 31, ^, |3q| . But in BIG around flat background, the tachyonic ghost cannot 
be avoided easily. There have been attempts to evade the problem by modifying the 
setup [13, |27], p3, but a satisfactory resolution is still lacking. 



3. The third issue is the mode of realization of the gauge symmetry, in this case 4- 
dimensional general covariance on the brane. This is related to the structure of the 
effective theory on the brane. It is known that in both Fierz-Pauli massive gravity, as 



well as in Higgs gravity setups (see [35, ^6[), a graviton mass is related to the breaking 
of 4-d general covariance. However, in BIG, the starting theory is manifestly invariant 
under bulk and brane gauge transformations, while at the same time 4-d gravitons 
acquire an effective mass. In usual treatments this result is obtained after gauge 
fixing which obscures the nature of symmetries of the effective theory. How is the 
generation of mass related to the realization of symmetry in BIG models? 



^Progress in this direction will be reported elsewhere. 



The fourth issue is a more detailed understanding of the technical aspects of the 
resulting modified gravity on the brane and its response to the cosmological constant, 
both in flat and curved background approximations. Technically, a zero thickness 



brane should consistently be replaced by one with an effective width [g, 17, p3| , 
before a quantitative investigation of graviton mass, its decay and other possible 
modifications of the gravitational force. It is also important to see explicitly how a 
cosmological constant affects the situation. In this paper such an analysis is carried 
out in the flat background approximation. 

5. Related to this is the fifth issue of naturalness of scales in BIG, at least based on 
some elementary considerations. The basic BIG model contains bulk and brane 
Planck scales, say, M* and Mp, respectively, as free parameters. For a realistic 
Mp ~ 10 19 GeV, to insure that the resulting IR modifications of gravity are not too 
large, the bulk gravity scale M* should be unnaturally low (< 10 _3 eV for n > 2), 
corresponding to strongly coupled bulk gravity pf. The effective brane thickness 



introduces another scale, though often this is taken to be related to M* [17, 23]. 
But a non-zero A also affects the identification of the gravitational constant and 
other standard model couplings, influencing the problem of scales. Taking this into 
account is important since as a viable theory BIG should remain calculable at the 
interacting level and should not generate new hierarchy problems. It is not obvious 
that BIG could meet these requirements. 

Considering the potential ramifications of brane induced gravity for the cosmological 
constant problem, we revisit these models and analyze them in some depth in the hope of 
gaining a better understanding of the issues outlined above. This work address some of 
these issues directly, and hopefully sheds light on, and sharpens the context of, the rest. 

1.2 Overview and discussion of our results 

This subsection contains a description of our results and their implications. We concentrate 
mainly on codimension n > 2 BIG models relevant to the cosmological constant problem 
and allow for the addition of extrinsic curvature terms and a brane cosmological constant 
(brane tension) A. In this model, we study linearized metric fluctuations around a back- 
ground configuration which is approximated as a flat 3-brane in a flat bulk spacetime (for 
the reason explained under item (1) of the previous subsection). The crucial point is that 
the presence of a small brane cosmological constant does not invalidate the flat background 
approximation. 

Spontaneous breaking of general covariance: The basic variables are the bulk 
metric Gmn an d the brane embedding functions x (<r). The symmetries are the brane 
and bulk general coordinate transformations, involving a^ and X M , respectively. The 
standard gauge choice (Monge or static gauge for the brane, and harmonic gauge for the 
bulk symmetries) is natural for bulk physics but obscures the nature of symmetries of the 
effective brane theory. Thus, with a focus on 4-dimensional physics, we work with gauge 
invariant variables and the induced metric g^ = GMNd^,x M d v x N , that naturally couples 



to brane matter. The outcome is that, beside the usual gauge invariant parts of the metric, 
we encounter 4 new gauge invariant variables. These are closely related to Stiickelberg 
fields used to restore gauge invariance in Fierz-Pauli massive gravity. The 4-dimensional 
effective action, obtained after integrating out all extra dimensional modes, depends on 
these variables, and through them, on gauge dependent components of the metric. This 
explicitly demonstrates how 4-dimensional general covariance is spontaneously broken by 
the bulk-brane setup in the effective 4-dimensional theory, whereas all symmetries are 
manifest in the starting action. Thus, BIG naturally implement a gravitational Higgs 



mechanism of the sort discussed in [|35|, |36[ . These issues are discussed in sections 6 and 7. 
Explicit analysis of graviton mass and decay: To study graviton propagation 
on the brane, one has to consider thick branes to avoid divergences associated with van- 
ishing brane width. It is then known that integrating out bulk modes contributes a "self- 
energy" term to the 4-dimensional graviton propagator, giving, schematically (see equations 

CTCTCT ), 

n (0,2) 

Bk 2 + Aoj n - 2 N(WE{k) ' ( ' ' 

A and B are inverse Newton constants in d and 4 dimensions, u> is the brane width and 
the superscripts refer to spin and 2 modes of graviton. The scale of graviton mass is set 
by Aui n ~ 2 /B, taken to be very small for phenomenological reasons. 

1. After reviewing the well known problem of the tachyonic ghost nature of the spin 
mode, we use a Gaussian profile for the brane thickness to explicitly extract from 
E(fc) a mass term, a decay width and a momentum dependent modification of the 
Newton constant for the healthy spin 2 mode. All this hinges on the fact that, as 
we show, the standard perturbative treatment of unstable particles in quantum field 
theory remains applicable for small graviton masses. 

2. The suppression of the probability for graviton escape into the bulk, resulting in 
effectively 4-dimensional gravity on the brane, is demonstrated using the optical 
theorem. 

Suppression of cosmological constant: In Einstein-Hilbert gravity, a cosmological 
constant A, however small, destabilizes flat space, giving rise to de Sitter or anti de Sitter 
spacetimes. In BIG, a brane A can curve directions transverse to the brane, leaving the 
brane metric flat. However, in the analysis of BIG around a flat background, commonly 
employed to extract more detailed information about brane gravity, the setup does not 
allow for absorbing A in the curvature of extra dimensions. But, as discussed in section 4, 
the theory compensates for this and neutralizes A using a S(0) ^ piece of the self-energy 
in the scalar mode propagator, closely related to its tachyonic mass. This leaves the brane 
metric flat 4 , and only modifies gravitational and non-gravitational couplings. Based on 
this analysis, the following observations should be made: 



4 The n > 2 BIG only attempts to explain why one does not observe a large cosmological constant. At 
least in this form, it does not explain the observed recent accelerated expansion due to a very small A. 



1. As far as filtering out A is concerned, the tachyonic mass of the scalar ghost replaces 
the function of a background curvature (as an artifact of flat background approxima- 
tion) and is not unhealthy. In fact, the mass is tachyonic only with respect to the 
wrong sign of the ghost kinetic term, but has a positive contribution to the energy. 

2. The momentum dependence of S(/c), which determines the graviton decay width and 
normalizes the gravitational coupling, has no bearing on the suppression of A (which 
is entirely due to S(0)). This is to be contrasted with the more general degravitation 



mechanisms [37] where a momentum dependent gravitational coupling could make 
the theory less sensitive to A. But the outcome of n > 2 BIG is closer to the scenario 
in |3ql, more closely modeled around massive gravity. 



3. In some of the literature on BIG, the suppression of A is heuristically explained in 
terms of 3-brane gravity becoming 4 + n-dimensional, hence weaker, at very large 
distances. This truly is the behaviour of n = 1 theory which, however, does not 
filter out A. A closer inspection of the propagator in n > 2 models shows that brane 
gravity at large distances remains 4-dimensional and is dominated by S(0), related 
to graviton mass, which suppresses A. 

4. A lesson from this analysis is that the flat background approximation to BIG remains 
valid in the presence of a small cosmological constant, with implications for the 
resolution of the ghost problem to be discussed below. 

Extrinsic curvature effects: Derivations of the BIG action in string theory setups 
also produce terms dependent on the extrinsic curvature of the brane. One may wonder 
if neglecting these terms has an effect on the ghost problem. The analysis in section 5 
shows that these terms have no effect on brane gravity at the linearized level. This may be 
disappointing from the point of view of the ghost problem but also shows that these terms 
do not drastically modify gravity on the brane. 

Status of the ghost problem: In equation ( |1.1[ ) for the graviton propagators in 
BIG, 0\j,i, pa contain the tensor structure. This is exactly the same as the tensor structure 
of Einstein-Hilbert gravity in 4 dimensions (obtainable by setting A = 0) . It is known that 
with this tensor structure, the theory is ghost free only for zero masses. Any non-zero mass 
pole implies a ghost, as is the case with BIG 5 . Resolving this problem is crucial for any 
eventual application of BIG. We do not achieve this in this paper, but sharpen the context 
of the problem. 

1. As a warm up consider the following puzzle: The BIG action can arise within seem- 



ingly consistent setups, like the string theory setup of [18]. How can these consistent 
setups lead to a ghost inconsistency? To answer this, note that the identification of 
ghost in BIG hinges on the S term in ( |1.1| ) being a small self-energy correction. This 



5 The ghost can also be avoided by changing the tensor structure gravity in BIG for example, as in |27], p4] 
but that introduces vDVZ type discontinuity in the free theory. These could be avoided in the non-linear 
regime by the Vainshtein mechanism [139, kd], but we will not consider such alternatives here. 



is obviously the case in the phenomenologically interesting regime of small graviton 
mass (~ Auo n ~ 2 /B) which requires that the brane Einstein-Hilbert action dominates 
over the bulk one. If this is not the case, then either the ghost may have a large 
mass, beyond the regime of validity of the low-energy effective action, or S(/c) may 
contain appreciable corrections to Bk 2 making the above interpretation of the prop- 
agator meaningless. This is the case in the string theory regime of [18| where, for a 



small string coupling g s , A ~ l/g 2 dominates over B ~ l/g s (assuming ui ~ g®), the 
propagator is essentially S _1 (/c), and gravity is ci-dimensional and ghost free. On the 
other hand, assuming w 4 ~ g s and n = 4 (more natural for D3-branes), the mass is 
at string scale. 

2. Could the appearance of ghost (in the phenomenologically interesting cases) be re- 
lated to the use of flat background approximation, ignoring A and its back reaction 
on bulk geometry? As discussed in section 4, since our analysis is already valid for a 
small A, one can answer this in the negative by a continuity argument: that in any 
case, the small A limit should correspond to our result which is not ghost free. 

One may formulate a sharper criterion: In a gravitational background sourced by 
a brane A, the metric along the brane remains flat, hence the background mostly 
affects gravitational dynamics in the bulk. In terms of the propagators in (1.1), the 



curved background mostly affects the ^4-dependent self-energy term E, leaving Bk 2 
and tensor structure unchanged. Then in the phenomenologically interesting regime 
of small mass, where the B term dominates, the propagator has the same structure 
as in flat background. The discussion above then shows that the ghost mode can be 
avoided only if the new self-energy term satisfies £ cu ™ed(0) = 0, implying no hard 
masses. This can be regarded as the constraint on admissible ghost free backgrounds, 
if any. For more on this see the concluding section. 

3. In the 4-dimensional effective action for BIG derived in section 7, the ghost can 
be related to the Stiickelberg-like fields that arise in the gauge invariant formalism, 
and hence to the spontaneous breakdown of general covariance. This underlines the 
similarity with the ghost in standard massive gravity and in Higgs gravity. 

1.3 Organization of the paper 

The rest of the paper is organized as follows: In section 2, we introduce the model and 
detail the setup for the analysis. Gauge invariant variables are introduced and a brane 
thickness profile is used to avoid divergences associated with zero-thickness branes. 

In section 3, we solve the equations for the basic BIG model and review the tachyonic 
ghost problem. The suppression of graviton escape into bulk is described in terms of 
the optical theorem. We then carry out an explicit analysis of the graviton propagator, 
extracting from it a mass, a decay width and a modulation of the gravitational coupling. 
It is shown in the process that standard perturbative QFT methods remain applicable for 
small graviton masses. 

In section 4, we consider the effect of a small brane cosmological constant A in the flat 
background approximation, showing that for codimension n > 2, it is filtered out by the 



tachyonic scalar mass. The tachyon problem is dismissed. The ghost problem is argued to 
survive, even when the theory is analyzed in a background sourced by A, except for special 
backgrounds (if they exist at all) that do not make brane gravity massive. 

In section 5, we solve the equations in the presence of extrinsic curvature terms that 
generically arise in derivations of the BIG action. It is shown that they have no effect 
whatsoever on the brane gravity. 

In section 6, we revisit massive Fierz-Pauli gravity with Stiickelberg fields in terms 
of gauge invariant variables. In particular, it is shown that the correct solutions can be 
obtained from a 1-parameter family of actions. 

In section 7, we derive an effective action for the brane gravity by integrating out 
bulk related modes. We construct a 1-parameter family of effective actions that reproduce 
the correct solution and, generically, contain Stiickelberg fields. This shows that gauge 
dependent modes of metric do not decouple, implying a spontaneous breakdown of 4- 
dimensional general covariance. The relation to the gravitational Higgs mechanism is also 
pointed out. The conclusions are summarized in section 8, and the appendix contains a 
brief discussion of gauge invariant variables and related zero mode issues. 

2. Preliminaries 

In this section we introduce the brane induced gravity model and describe our setup for 
its analysis: the use of gauge invariant variables, description of thick branes in terms of a 
density profile and its use as a regulator. Finally we consider the bulk-brane relations and 
the appearance of 4 new gauge invariant variables peculiar to the bulk-brane setup. 

2.1 The model and its origins 

In brane induced gravity models, the Universe is regarded as a 3-brane with coordinates 
c M (fj, = 0, • • • , 3) embedded in a non-compact (d = 4 + ra)-dimensional bulk spacetime with 
coordinates x M through x M (<r). In the basic model, dynamics of gravity is governed by an 
Einstein-Hilbert action Sgjf for the bulk metric Gmn{x) and a similar action Sgfp e for 
the brane induced metric, g^ v = d fl x M d u x N GMN, sourced by brane matter S^ ane |l], [2|. 
In this paper, we also allow for adding a brane tension (brane cosmological constant) term 
gbrane an( ^ ex t r insic curvature corrections Sq ane , to be specified later. Then 6 , 




GRW -B f d 4 a^g~R^ + S% ane + S% ane + S b ™ ne . (2.1) 



As argued in Q ||, 17 1, such models naturally arise in braneworld setups as a result of 
integrating out massive matter on the brane following [H3]. More specifically, they were 
shown to arise in non-supersymmetric string theory, as a result of closed string scattering 



off D-branes fig , 19 1 where the action fl2.1| ), including S^ ane can be computed explicitly . 



6 We follow the sign conventions of Weinberg in mU 

7 The calculation of [TLSl, llSl] leads to the expected behaviour of A ~ g~ 2 and B ~ g^ 1 with the string 
coupling g s . Mass dimensions are supplied by the string tension, 1/a' . A curious feature of their result is 
that while, A > 0, one gets B < 0. Also, the relative sign between S^x™ e and Sf[ ane is opposite to what 
appears in the Gauss-Codazzi equation. However, for the BIG results to be relevant to phenomenology, we 
must take B > and adjust the parameter values appropriately, as will be specified later. 



At this order in the ol perturbation theory, there are also bulk R 2 terms, in the Gauss- 
Bonnet combination, that do not contribute to our analysis, as well as terms involving 
other massless string states that are ignored here. A different string theory realization was 
proposed in [20]. Here for the couplings A, B, etc. we do not use their calculated high 



energy values, but treat them as low energy phenomenological parameters. 

The d = 5 model is the well known DGP model that has been extensively investigated. 
As argued in []|, |6| it cannot address the problem of the observed smallness of cosmological 



constant, d = 6 is also special and was investigated in ||24| , |25| , |26|| . Below, we are interested 
mainly in d > 6 models. 

2.2 Flat background and gauge invariant variables 

We investigate brane induced gravity at the linearized level around a flat background, 
which is technically the easiest and a good first approximation. Consider a flat bulk with 
a flat brane and split {x M } into {xft ,x l ± }, respectively, parallel and perpendicular to the 
brane. The flat brane corresponds to xft = a 11 , x l ± = y l Q (constants). Linearizing Gmn(%), 
9iiv{&) and x (a) around this background gives the fluctuations HMN,h^ u , f 1 and y l , 

Gmn =Vmn + Hmn(x) , g^u = VtJ.u + h flu (a) , (2.2) 

xfj(cr) = a» + /» , x\{a) = yj + y\a) . (2.3) 

To first order, h(x(cr)) = h(xn), etc., and the fluctuations are related by the linearized 
pullback equation h fJLU {x\\) = -ff MJ ,(x||,yo) + dpfv + d u f^. Infinitesimal bulk and brane 
diffeomorphisms transform the fluctuations as, 

6x m = ^m . 6 ^ Hmn = _2d (M £ 7V) , Szhp, = o , stf* = e , 

5a» = \» : 8 X H MN = , <J A V = ~ 2d ^K) , S\P = -A" . 

The last equation arises since x M (a) are scalars under brane diffeomorphism so that 5\f^ = 
Sxx^ = —\ v 5y. The ffj, play a role similar to Stiickelberg fields in massive gravity. 

To solve the equations of motion, we use gauge invariant variables instead of gauge 
fixing transformations ( p.4[ ) 8 . This disentangles the spin and spin 2 modes of the graviton 
and promotes /„ to new gauge invariant variables. The bulk field Hmn decomposes as 

H M n = H MN + d M A N + d N A M + d M d N $ + - t] M n S (2.5) 

d 

where d M H MN = 0, H M M = 0, and d M Am = 0. Then, H MN and S are gauge invariant 



while Am and $ give rise to the transformations of Hmn m (kM)- The relevant projection 
operators are listed in appendix A. The bulk Einstein-Hilbert action becomes, 

S^ k = - j f d d x(d L H MN d L H MN - d L Hd L H + 2d L H8 N H LN - 2d M H MN 8 L H LN 



-j i ,r.v 



,L u ±MN a „± (d-2)(d -l) 



d L H ±MN d L H MN - ± ^ J -d L sd L s 



(2.6) 



3 Customarily, one fixes static (Monge) gauge, / M = 0, on the x^(a) and harmonic gauge on Hf. 



where gauge dependent terms drop out. To obtain the contribution to the equation of 
motion from the second line, express the variations 5H^ M and SS in terms of 5Hmn using 
the projection operators in Appendix A. The linearized Einstein-Hilbert operator becomes, 

(£dH)MN = y^dH MN H — — — (8m9n — VMN^d) S . (2.7) 

Similarly, the induced brane field h^ is decomposed as 

hpv = hj; v + d^a v + dydy, + d^d v cj) + - r?^ s , (2.8) 

where d^h^ = 0, h^ = and d^a^ = 0. hh v and s are gauge invariant. The brane 
Einstein-Hilbert action 5^g" e and operator (£±h)n V can be read off from ([O]) and ( EO| ) for 
d = 4, after appropriate field replacements. There are also new gauge invariant variables 



based on ft* and y l to be introduced later in equations 2.17 and f>. 7j . 

2.3 Thick branes and "blurred" quantities 

Solving the equations of motion involves the massless scalar propagator in d dimensions, 

„ ,. ik(x\\— x',.)+iq(xj_— x',) 

G(x ll -x\ l ,X ± -X l ± ) = -J d'kj d V k2+q2 _ ie • (2-9) 

kfj, and qi denote momenta parallel and transverse to the brane. For two points restricted to 
lie on a zero thickness brane, x\ = x']_ = y l and the q integral in G(xn — x',,,0) diverges for 
n > 1. Correspondingly, a bulk field S sourced by a brane localized source 5(x± — yo)T(x\\), 
i.e., S(x\\,x± — yo) = J d x '\\G(x\\ — x',,,x± — yo)T(x\\) also diverges as x± — > yo- Hence, 
the restriction of such bulk fields to the brane is ill defined. Essentially, a zero thickness 
brane is a point source in transverse directions leading to a \x± — yo|~ divergence, as 

in the Coulomb or Yukawa potentials. 

This problem is resolved by realizing that dynamics give branes a form factor and 
hence an effective thickness, as discussed in [p3| l for D-branes and in |17] for solitonic branes. 
This is taken into account by replacing S(x± — yo) by a brane thickness profile function, 
P(x±—yo), with normalization J d n xj_P = 1. Then, the divergent brane restriction S(x\\, 0) 
of the bulk field is replaced by the "blurred" field (S , )(xm), a weighted average over the 
brane width, 

(5)(x,|) = J d n x ± P(x ± - yo)S(x\\ t x ± - y ) , (2.10) 

With Fourier transform conventions P(x± — yo) = f d n q P(q)e iq ^ Xj -~ yo ', one gets the mo- 
mentum space equation 

(S)(k) = (27r) n Jd n q P(q)S(k,q). (2.11) 

In the same way, the divergent brane restricted propagator G(x\\ — Xm,0) gets replaced 
by its well defined thick brane analogue, 



{G)(x\\-x'u)= I d n x±d n x' ± P(x±-y )G(x\\-x'u,x±-x' ± )P(x' ± -y ) 
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Physically, this is the propagation amplitude from xn to x'u located on the brane, with 
a small derealization in the transverse directions. This does not yet take into account 
physical effects due to the presence of the brane. However, this construct will contribute 
to the physical propagators and the interesting physics associated with it will be discussed 
in the next section. 

The corresponding momentum space expression (with P depending only on q 2 ) is 

(g)(AQ = -(27r)" J 'd n g k } P + { f_ ie - (2-12) 

Again, for a zero-thickness brane, P = 1 and the expression diverges for n > 1. For 
thick branes, P(q) effectively implements a UV regularization at high q through a length 
scale associated with the brane width, u. Earlier, a non-zero brane width has been used 



to justify implementing a sharp cutoff ~ 1/u on the q integral [23, 27], or solving the 
equations separately outside and inside the thick brane [O, |34|. Here, we will explicitly 
retain a smooth profile function P through which the results depend on the brane width 
uj. Such dependences are not affected by the actual form of P, that encodes high energy 
effects and which will only change the numerics. 

Since the normalized P has mass dimension n, P(q) is dimensionless. Then rotational 
invariance in the transverse space means it only depends on the combination q 2 uj 2 . This 
implies (after simple manipulations in ( 2.12; )), 



(G)(k) = -^Z-\u; 2 k 2 ), (2.13) 

depending on k only through u = uj 2 k 2 . The expression captures the form of divergence as 
uj — ¥ 0. Later, S will appear in the same way as a self-energy correction in quantum field 
theory, hence the notation. Note that w n {G) has the correct dimension for a propagator. 

The use of the brane effective width uj means that the theory is valid for describing 
interactions of transverse gravitons with the brane as long as q < 1/uj. Probing the brane 
at shorter scales is meaningless in the low energy theory. To insure that interactions 
too cannot probe beyond the brane width, the effective theory must also be restricted to 
u = (kuj) 2 < 1, for momenta along the brane. 

For the most part we do not need the functional form of P. But for explicit calculations, 
a natural choice, motivated by [2C], is the Gaussian form 9 , 

P(x±) = ^ e -(^/ 2 -) 2 , P(q) = -^— e -i 2 " 2 n . (2.14) 

The profile function can also be discussed in a covariant setup, but that is not needed here. 

2.4 Bulk-brane relations 

For a thick brane, the pullback equation relating bulk and brane metric fluctuations be- 
comes, 

V = <*W + d nf» + ft//n . ( 2 - 15 ) 



9 In contrast with a sharp cut-off on q, a Gaussian P(q) leads to a non-negative P(x± — j/o), consistent 
with its interpretation as a thickness profile. The main results are not affected by these choices. 
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with only x\\ dependences. Note that fields intrinsically defined on the brane are not 
affected by the "blurring" procedure ( 2.10J ). Then, using (|2.5| ) and ( |2.8| ) gives, 



(H L ) 



fill 



'III' 



d.F u - a v F, 



'///;■ 1 j(S) 



(2.16) 



The F^ are new variables invariant under both bulk and brane gauge transformations, 



f fi + (A ti )-a ll + -d^) 



(2.17) 



The gauge dependent variables in Fn do not appear in the action (2.1). However, solving 
the equations of motion requires either gauge fixing or using gauge invariant variables along 
with ( |2.16| ). Thus, in bulk-brane setups, gauge variant variables survive through F^ and 
contribute to the solutions for the brane fields s and h 1 - . The implication of this as a broken 
phase realization of 4-dimensional general covariance is discussed in the last section. 



3. The Basic Brane Induced Gravity Model 

In this section, we review and further investigate the solutions of the basic BIG model 
HI, H based on bulk and brane Einstein-Hilbert actions of ( |2~T| ) but without £* rane and 
Sa ane . First, we solve the linearized equations of motion for the brane fields and review 



the tachyon/ghost problem [28|. The propagator is then analyzed in depth, extracting a 
mass, a decay width and a modulation of the gravitational coupling. The suppression of 
graviton decay into the bulk is explained with the help of the optical theorem. The reader 
not interested in the details of solving the equations of motion can directly jump to the 
solutions (|3.15| ) and ( |3.16| ). 



3.1 Equation of motion and solutions 

For the basic BIG model, the linearized equation of motion for Hmn expressed in terms 



of Einstein-Hilbert operators (|2.7| ) and for thick branes, takes the form 



A(£ d H) MN + P(x x - y ) ( B(£*hr + -T"" ) <<5 



1. 



M X N 



0. 



(3.1) 



In our momentum space conventions, it becomes 10 (with p M = {k^,q 1 }), 

d-2 



-A 



{k 2 + q z )H ±Mn + 



d 



( v M v N_ r] MNnJ2 



(k 2 + q 2 )) S 



-p(q) 



B I k 2 h^ u + - {k^k u - rj^k 2 ) s 



(fc,«) 

rpflU 



W 



6 M 6 N 



0. 



(3.2) 



This is to be combined with the surface equation ( p. 16 ). In the following, the ie terms 
in the bulk propagator are not always written explicitly, but are finally included in the 
blurred propagator ( p. 12 ). 



3 The (fits) components alone can determine the brane fields, but it is more convenient to use all equations. 
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The timn -trace of (|3.2|) determines S in terms of s, 

d P(q) 



S(k,q) 



A(d-l)(d-2) k 2 + q 2 -ie V " 2 



Tf! + lBk 2 s 



(3.3) 



(fc) 



Restricting to the brane by using (2.11) and the expression for (G) in ( |2.12| ), one gets, 

d ,-=, /-„ 3 



<S>(*) 



(G) (r; + ^fe 2 s 



A(d-l)(d-2) x ~' V " ' 2 
For later convenience, we express S(k,q) in terms of (S)(k) 

P(q) 1 



W 



(3.4) 



S(fc,g) 



k 2 +q 2 (G) 



(S). 



(3.5) 



T/ie (p, i) components give (since "blurring" now involves integrating over an odd 
function of q l ), 



H^ l {k,< 



d-2 k^q i 



S(k,q) 



{H L )^{k) = 0. 



d k 2 + q 2 

However, note that using ( |3.5[) , and the transversality of H , one gets, 

d-2 



{ qi H^)(k) = -k v (H- 



\fW 



■k^iS) 



d (G) 

where we introduce a new function (g)(k) (that will drop out of most expressions), 

q 2 P 2 



(3.6) 



(3.7) 



(27T)" J 



d" 



(k 2 



.,2^2 " 



The (i,j) components give, 
H ±ij (k,q) 



d-2 , q % qi 



V iJ )S(k,q). 



(3.8) 



(3.9) 



d k 2 + q 2 
From this one can immediately see that (qiH J ) = and, for i ^ j, (H • J ') = 0. However, 



<#/ J }(*0 



*Z* ® + „ (S) , < ra ^«) W = -<Z» * ® 0) . (3.10) 



The surface equation ( p. 16 ), on taking a trace and a double divergence, gives 
4 



(H^) = -2ik»F^ ~^(S) + s, k^k"(H L )^ = -k 2 UiWFn + i(5> - J?) . (3.11) 

Since (H^) = -(H^) and k^k v {Ji^)^ = {q l q j H ± ^){k), on using flJUOD , (|37TT|) become, 



2JFF,, 



(d-2)(-^-+n)-4 

(G) 



a 



(d-2)-^--l 



f + I. (3,2) 



13 



Eliminating k^F^ gives a relation between (S) and s as, 

(d-lKd-5) (§)= _3, (313) 

Substituting back in the expression for k^F^ determines it in terms of s. Combining this 
with the divergence of ( |2.16| ) and then using (|3.7| ) leads to the solution for the F^ fields, 

iW „ 1 d-2 ( ., , 

"" = -- iT ,Vs, where, V = - _ _ In - 3 -^ , CM,, 




Finally, combining ( |3.13 ) with ( |3.4D gives the solution ( [3.15 ) below for s in terms of T = Tfc. 



The (ji,v) Components of Q3.2| ) can be solved for H^ v . One can then compute (H^ v ) 
using ( |3.5| ), Using ( p. 16 ), this becomes an equation for h~; u in terms of S, J and F^ all of 
which are known in terms of T^ v . Thus one has the final solutions, 

^ = -in k 2 + A L ^ r, (3-15) 

V = s fea.A^-! (^ " 3 { ^ -^T )T )- (3 - 16) 

From this one can directly read off the gauge independent parts of the Greens functions 
G ' , v ,, or the 4-dimensional gauge invariant amplitude, 

d A x I d 4 x'T» u (x)G^ v ,(x - x')T^ u '{x') 

= J d 4 xT^(x)h^(x) = j'd 4 x (T^h^ u + \TsY (3.17) 

If needed, the h^ in a given gauge can be constructed by adding gauge transformations, 
h/j,u = h^ u + \ ri^ s + d^v + d^, and solving the gauge conditions for the £ M . 

One can immediately draw a number of conclusions, mostly known in the literature, 
based on the structure of s and h 1 - solutions: 

(1) Zero thickness limit: Note that in the absence of extra dimensions, (G) n =o = 
— l//c 2 , leading to the solutions for ordinary Einstein-Hilbert gravity, 

(^)(n=o) = -j, p (t^ - -(v^u ~ -p")rj , (i)n=o = - ^7 p r • (3.18) 



Also (gr)n=o = gives (i^) n=0 = 0. The BIG results (gig), (gig) differ from this only in 



the (G) _1 dependent terms. Now, while the brane width u was introduced to keep (G) 
finite for n > 1, the final solution is well defined for uj — > 0. Then, (G) diverges (from 
( 2.13 )), resulting again in ( 3.1&J ). Hence, for n > 2, modifications of 4-dimensional gravity 



arise only for thick branes. Since the spin 2 and spin parts of the propagator in BIG 
have the same coefficients as in ordinary gravity (except d = 5 where s = 0), there is no 



vDVZ discontinuity |3lL 32, p3[ in the u — > limit 
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(2) Tachyon problem: For n > 2, (G)(k 2 = 0) is finite and contributes a hard 
mass to the propagators. The scale of the mass is set by Auj n ~ 2 / B and the gravitational 
coupling is given by Gn ~ 1/B. Thus, for phenomenological reasons, B should be large 
and Auj n ~ 2 should be small. But (G) _1 appears with opposite signs in the denominators, 
indicating that if h has a healthy mass, then s will be tachyonic or vice versa jp8| . For 
the conventional choice of B > and A > 0, h cannot have a tachyonic pole (since for 
any k 2 > 0, ( p,12| ) gives (G)(k) < 0). Then, s is tachyonic. 

3) Ghost problem: It is well known [pl| [44| , |45|1 that a propagator for h^ u with the 
tensor structure of ordinary massless gravity in 4-dimensions, is ghost free only for zero 
graviton masses. Thus the non-zero masses and the —2/3 factor in the s solution makes 
the tachyonic spin mode also a ghost J28|] . The presence of this tachyonic ghost in brane 
induced gravity has so far hampered its further development. One of our purposes in this 
paper is to shed light on the origin of this ghost mode. 

4) Unstable gravitons: Ignoring the tachyonic ghost s field, the massive spin 2 
graviton h 1 - on the brane is not a stable particle but a resonance state. Technically, this 
is because (G) has a branch cut from k 2 = to — oo due to the continuum of Kaluza- 
Klein modes resulting from an uncompactified bulk (|j, g]). For very small masses, the life 
time is very large. In the next subsection we will explain the origin of the decay and its 
suppression using standard quantum field theory and then carry out a detailed analysis of 
the propagator. 



5) String theory limit: The string theory computation of the BIG action in [[j^, 19] 
obtained B < 0, making h 1 - the tachyonic ghost. But in this setup, in terms of the string 
coupling g s , A ~ 1/g 2 and B ~ l/g s - Assuming uj ~ g®, weak string coupling corresponds 
to large masses, beyond the validity limits of the theory. Even for D3-branes where n = 6 
and w 4 ~ g s is a more natural choice, masses are ^-independent and at string scale. In 
either case, for k 2 values within the validity range of the theory, the propagator is given 
by —(G)(k)/A, avoiding the tachyon/ghost problems, and gravity is essentially higher 
dimensional. However, the orbifold based construction of p(| that attempts to derive BIG 
from string theory with the phenomenologically interesting parameter ranges, does not 
evade the tachyon/ghost problem in this manner. 

In the rest of this section we concentrate on the healthy spin 2 mode, returning to the 
origins of the ghost mode in the next section. 

3.2 Graviton decay from optical theorem 

A main feature of BIG, analyzed around flat background, is that h is a massive unstable 
graviton in 4 dimensions due to the possibility of gravitons escaping into the bulk j|, S, 0. 
However, the graviton lifetime must be large for phenomenological reasons. Below, we 
show that the suppression of graviton decay can be explained using general field theory 
arguments. The detailed structure of the propagator is analyzed in the next subsection. 

The bulk propagator G(x — x') (£j|) gives the amplitude for graviton propagation 
from any x to any x' (ignoring the tensor structure and for canonically normalized fields, 
so no 1/A factor). But (G)(x\\ — x',,) is constructed to describe a restricted propagation 
between two points with their x± coordinates delocalized over the width w of a thin region 
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that would be occupied by the brane. Of course, there is a finite probability a escape that 
gravitons emitted at x\\ within the thin region, end up somewhere in the bulk. Such 
processes are not described by the restricted propagator (G) and appear to it as decay 
channels, thus giving it an imaginary part. The relationship is quantified by the optical 
theorem where, in terms of the corresponding amputated 2-point function to~ n (G)^ 1 (with 
u factors restoring the correct dimensions for the propagator), 

2Im(u- n (G)- 1 )~a escape ^0. 

This is a geometric result that holds for any region of space that would be occupied by a 
thick brane and applies to gravitons not yet dynamically affected by a physical brane. 

Now, with a physical brane in place and after taking gravitational dynamics on the 
brane into account, we obtained a brane-to-brane propagator Gbb = — {k 2 — ■r(G) -1 ) -1 



in ( 3.16 ), but now written for the canonically normalized brane field and with suppressed 
tensor structure. The stability of the corresponding state is again encoded in the imaginary 
part of the amputated 2-point function which now gives, 

~ A ~ Auf 1 

21m G^ 1 = 2- Im {G)' 1 ~ — a escape . (3.19) 

Hence, graviton decay on the brane is still due to escape into the bulk, but with an ampli- 
tude suppressed by a factor Aco n /B as compared to the unhindered escape in the absence 
of brane gravitational dynamics. This is a way of seeing that a large brane Einstein-Hilbert 
term confines gravitons to the brane by suppressing their escape into the bulk. 
Finally, for completeness, one can also compute the brane-to-bulk propagator, 



S»(M)-- PW(G >" 



k 2 + q 2 



1 1 



B k 2_A {G) - 

The brane-to-bulk decay probability (Tbrane^buik can be constructed out of this after proper 
amputations (using Gbb for the brane leg and G for the bulk leg). One can then verify the 
optical theorem for the physical propagator, 

21 m G hb ~ Obrane^rbulk 

In this sense, the verification of the optical theorem on the brane requires using the brane- 
to-bulk propagator, since the brane theory alone is not unitary. 

3.3 Graviton mass, decay width and coupling 

Let us investigate the properties of the spin 2 graviton h?~ ( |3.16| ) in more detail. The physics 
of the corresponding brane-to-brane propagator follows from the structure of (G)(k 2 — ie). 
We will see that the propagator has the standard form for a massive unstable particle in 
QFT. More importantly, for small masses it is amenable to the standard particle physics 
approximation methods. 

To evaluate (G) ( p. 12 ) we use the Gaussian form ( 2.14 ) for P(q). Doing the angular 



integrals in (/-space gives, 



n _! < ■ ..-_, . »_,n. 



((V) - -2.V / dqq n - i k2+g2 _ ie , N~ l = (4tt)* T(-) . (3.20) 
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First, let us consider the even n case and set 2m = n — 2 (the result for odd n is given 
later). After some manipulations, 



(G)(k) 



(-l) m N 



UJ 



n-1 



U 



3=1 \- u ^-< 



(3.21) 



where u = u k and we have used the notation [46|, 



E\(u — ie) 



e"' 



dt- 



(3.22) 



Since d u E\{u) = —e u /u, for £ 1 = uj n 2 {G) one gets 



E^iu-ie) 



-l) m N 



u m e u E 1 (u - ie) + $^(-l) r (r - 1)! 



u 



r=l 



(3.23) 



The corresponding expression for odd n is given by equation (3.35) below. 

Remember that the ie prescription tells us how to handle the k 2 + q 2 = poles in 



(2.12). Now, after performing the ^-integrations, the same prescription will dictate the 
correct e — > limit. To see what happens, note that the function Ei(z) is the analytic 
continuation, to the complex plane, of the exponential integral Ei(—u) defined, for real u, 
by 



Ei(-u) 



oo e -t 



dt- 



t 



(3.24) 



Ei(—u) is defined over the entire real line, but E\(z) has a branch cut from z = to 
z = — oo along the negative real axis. The crucial point is that as e — > 0, we approach the 
branch cut from below and 



lim£'i(n — ie) = —EA—u) + iir 9(— u) 



(3.25) 



where, 0(—u) is the unit step function. The imaginary part is the origin of the complex 
pole in the propagator, and hence of graviton decay . Writing S _1 = a + ib, the real and 
imaginary parts are given by 



a(u) = (-l) m N 



u m e u Ei(-u) - Y^(-l) r (r - 1) 



r=l 



b(u) = -7T iV \u\ m e" |u| 9{-u) 



(3.26) 

(3.27) 



Finally, in terms of £ = Si +iY,2, the brane-to-brane propagator in the h solution ( 3.16J ) 
is (with canonical field normalization and suppressed tensor structure), 



Gbb 



k 2 



B 



£i(w 2 A; 2 ) + iZ 2 {u 2 k 2 ) 



(3.: 



"Approaching the branch cut from above (corresponding to Ei(u + ie)), would give —in which does not 
lead to decay. So the correct sign is dictated by the ie prescription in the bulk propagator. 
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The above analysis then gives, 

£i(u) = ^2, S 2 (n) = -^>0. (3.29) 

This propagator has the familiar quantum field theory form where the analogue of £ 
arises due to self-energy corrections. It has a complex mass pole and, strictly speaking, the 
graviton mass and decay width are given by the real and imaginary parts of the complex 
pole. In practice, determining the location of the complex pole is not easy. But since the 
determination of mass anyway becomes somewhat arbitrary due to the finite life-time, an 
approximation can be made provided, near the mass pole, |£i| is much larger than 1 5^2 1 - 
In such cases, a mass can be defined as the pole position with only £1 present, while £2 
gives the decay width (see, for example, 47, |4^]). This is often a valid approximation in 



perturbative quantum field theory. Fortunately, one can show that it also holds here for 
small enough masses: For m > 1, as u — ¥ 0, £2 vanishes as u m while £1 goes to a constant 
— l/[N(rn — 1)!]. Hence, for a small enough mass, at the mass pole M 2 = —u/uj 2 , one has 
|£i| >> I £2 1, as wished. 

Then, in this approximation the mass M 2 = —k 2 is determined, on ignoring £2, by 



B 



u z k z = £i(urr). (3.30) 



Au r 

For n > 2, £i(0) - ~[N(m - l)!]" 1 < 0. So, by continuity, £i(w 2 /c 2 ) < for small 
arguments in the neighborhood of zero. Then the graviton is non-tachyonic, M 2 > 0. 
More concretely, this is visible in Fig. 1(a) where the right hand side of the above equation 



is plotted against u = <jj 2 k 2 for n = 6, 5, • • • , 3 {odd n curves are based on (|3.35|) ) and the left 
hand side is plotted for two different slopes B/Auj n , with values chosen only for illustrative 
purposes. The intersections determine the masses M 2 = —u/lo 2 , with larger slopes resulting 
in smaller masses. The plots also show that £1 is monotonic for n > 2 and the propagator 
has a single mass pole. In Fig. 1(b) we plot £2 vs co 2 k 2 . Comparing curves for the same 
n in both plots again shows that for very small masses, |£2( — lu 2 M 2 )\ << |£i(— co 2 M 2 )\, 
justifying the approximation (the curves for n = 1,2 look rather different and are not 
shown here). 

Now, in terms of the mass as determined above, the propagator has the standard form, 



Z(k) 

k 2 + M 2 - i 

where, 



G »- k 2 + M i-i^Z(k)Z 2 (k)' <3 ' 31) 



^..- ^ww; . (3 , 2) 

This will look more familiar when £i(w) is Taylor expanded around u = —u 2 M 2 . 

Note that Z retains a momentum dependence which amounts to a momentum depen- 
dence of the effective gravitational coupling, 

G e J f {k) = Z(k)/B. (3.33) 



18 




(a) Determination of Mass 



-150CH 

(b) Decay Width 



Figure 1: Behaviour of £i(lt) (left) and ^(w) (right) for codimensions n — 6 (the lowest curves) 
through n = 3 (close to the x-axis) 



This momentum dependence can be expressed entirely in terms of the dimensionless pa- 
rameters u and Au n /B. In particular, for time independent sources (k 2 > 0), S2 = and 
this is the only modification besides the graviton mass. 

A simple analysis of Z(k) shows that at large momenta, \k 2 \ » M 2 , Z ~ 1. But for 
small momenta, \k 2 \ ~ M 2 , on Taylor expanding around the mass pole given by ( |3.30| ), one 
gets Z < 1. Hence the gravitational coupling is weaker at large distances, but the variation 
is not drastic. This is over and above the effects due to graviton mass and decay width. 

Now, comparing ( 3.31| ) to the relativistic Breit-Wigner form near the mass pole, the 
decay width is identified as 



r = 



Au n - 2 



Mfi- Z(n)S2(n; 



= -MZ(u)^± 



>0, 



(3.34) 



u=-uj 2 M 2 



Once again, note that it is the ie prescription in the bulk propagator that leads to the 
correct sign for the decay width. 

For odd n the analysis proceeds along similar lines, except that the exponential integral 
is replaced by the error function and, with 2m' = n — 1, we get, 



TT\u-i€) = -(-l) m A% 



m'-i 

u 2 e 



u (l - erf(Vu^ie) 
4lHr^^l. (3-35) 



r=l 
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In fact, this expression has been used to 
plot the curves for n = 3,5 in Fig. 1(a) 
and Fig. 1(b). A unified description of both 
even and odd n is given in terms of the 
Whittaker functions, 



Ui 



'ea 



(4tt) 



W_n n^ (u — ie) 




0.002 



For completeness, in Fig. 2, we show a 
representative plot for the spectral density 
function, 

p(s) = ImG bb (k 2 = -s) , 

to show that it too has the standard form 
for an unstable particle, again with the 
Au n /B value chosen only for the purpose 
of illustration. 

To summarize, we have shown that although the origin of the "self-energy" contribu- 
tion, X, to the propagator in the BIG setup is very different from its origin in perturbative 
quantum field theory, the resulting unstable massive particles behave in very similar ways 
and are amenable to the same approximation. Hence, the modification of the spin 2 gravi- 
ton propagator ( 3.16J ), contained in (G) -1 , breaks up into 3 contributions: a hard mass M 
given by ( 3.3C ), a decay width T (|3.34|) and a momentum dependent Newton "constant" 



Figure 2: A Sample Spectral Density Function 



(|3.33[) . 



4. Screening of A and the Origin of Tachyon/Ghost Problem 

The possible resolution of the cosmological constant problem is the main reason for interest 
in n > 2 brane induced gravity, despite its ghost/tachyon issues. Support for this expec- 
tation comes, as discussed in || g], from the classical solutions of (i-dimensional gravity 
jlj, [f^] that describe a 3-brane with an ADM mass density corresponding to a non-zero 
brane tension (which is the same as the brane cosmological constant). To be precise, the 
classical solutions describe this setup at distances far away from the location of the brane 
core. The interpretation is that, for codimensions n > 2, a brane cosmological constant A 
can curve directions transverse to the brane rather than those parallel to it, keeping the 
brane metric flat 12 . Thus the effect of A is absorbed by bulk curvature and it does not 
curve the brane worldvolume to a de Sitter spacetime. 

Now, the presence of a brane Einstein-Hilbert term in the action does not affect this 
classical solution. But for a large enough coefficient, it is expected to dictate the dynamics 
of metric perturbations huv around the classical solution. This is why one expects n > 2 



12 14 p| a ^ so consider time dependent solutions but these do not show up in our linearized analysis and 
are not discussed here. 
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BIG models to describe effectively 4-dimensional gravity with suppressed A (The n 



model is somewhat special and has been analyzed in some detail in p4| , pq , 26]). 

Unfortunately, implementing the above scenario is not as straightforward. In reality, to 
study metric perturbations huv, one needs classical solutions that are valid in the vicinity 
of the brane (at least with some effective thickness included). The existing solution [14, 15] 
is not adequate for this purpose since the energy momentum tensor that sources it near 
the brane core, does not correspond to a brane A source. 

Hence, in practice, in BIG with codimension n > 2, the analysis of metric perturbations 
has mostly been performed in a flat background, instead of one sourced by A (even though 
the latter is used to argue for the suppression of A). It is also in the flat background 
analysis that one encounters the tachyonic ghost problem reviewed in the previous section. 
One may conjecture that this problem is an artifact of the flat approximation and would 
go away once a full fledged analysis is performed around a curved background sourced by 
A. However, in the absence of a classical solution adequate for such analysis, here we point 
out that the flat background analysis of the previous section is, by itself, sufficient to study 
the response of gravity to a small brane cosmological constant and exhibit its suppression. 
This also has ramifications for what to expect from a curved background analysis for the 
resolution of the ghost issue (as will be detailed below). The healthy spin 2 part of the 
graviton, on which we have concentrated so far, is not relevant here. Below we return to 
the unhealthy scalar mode to find that, 

• The "tachyonic" mass is actually healthy and is instrumental in filtering out A in 
the flat background approximation. It arises, at least partly, as an artifact of the flat 
approximation and as a substitute for the background curvature. 

• For phenomenologically interesting values of parameters (corresponding to a small 
graviton mass) , the ghost nature is more intractable and will most likely persist even 
when metric perturbations are analyzed in the background curved by A, unless one 
ends up with massless gravity on the brane in a fully non-linear background. 

• The 4-dimensional Newton constant and other couplings are modified by the screened 
A, consistent with y, ||[ ||. This may help in avoiding a very low bulk gravity scale 
A, but may introduce hierarchy problems for other couplings. 



• 



Finally we emphasize that BIG and massive Fierz-Pauli gravity (as studied in [pq ]) 
use different mechanisms to respond to A, in spite of a superficial similarity. 

4.1 Filtering out of brane cosmological constant 

It is easy to see that brane induced gravity for n > 2 is not very sensitive to a brane 
cosmological constant: The addition of the brane tension/cosmological constant term, 
5a = —A J d 4 a s/—g(X(a)), to the action amounts, at the linearized level, to the shift 

T^(k) -> f'^{k) = T^(k) + kr ] ^5^(k) . (4.1) 

The validity of the perturbative analysis around flat background requires that A is small, of 
the same order as T^. In ordinary Einstein-Hilbert gravity, A, however small, destabilizes 
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the flat space solution toward a de Sitter space (the manifestation of this at the propagator 
level is described below). Hence, in a modified theory of gravity the stability of flat space 
in the presence of A is an indication that gravity has been made less sensitive to the 
cosmological constant. The analogue of this argument for Fierz-Pauli massive gravity is 
considered in |jj|, although the actual mechanisms differ. 

Let us write the solutions ( 3.15J ) and ( |3.16| ) explicitly in terms scalar and traceless- 
transverse components of T^ v , 

2 1 ~ ~ 1 1 ~ 

m = -W k 2 + A_^ {d) -i Ts > h ^ = B & - A (G)-^ - (42) 

As discussed in Appendix A, it is only in this form that the solutions are valid for a 
cosmological constant source. For such a source, T^ = and T As = 4A<r 4 )(/c), leading 

to *^U = and, 

s A (k)=4c5^(k), where c = -| ^ ^ (G)(0) . (4.3) 

The crucial point is that only for n > 2, (G)(0) is finite and the analysis well defined 
(for example, for a Gaussian profile and for even n, (G)(0) = — (§ — 2)!/[(47r) 2r(|)w n_2 ]). 
Thus, in coordinate space, A effectively shifts the background from s = to a constant value 
s = 4c, keeping the flat background essentially unchanged. It is instructive to contrast this 
with Einstein-Hilbert gravity where one would get sa ~ A5^ 4 '(k)/k 2 , which is a solution 
to IZUsa ~ A. Then s A is quadratic in distance and diverges at large distances, indicating 
that flat space is not a good starting point for a perturbative expansion of de Sitter space. 
In contrast, in BIG, A does not produce a non-flat metric on the brane and is hence filtered 
out 13 . 

The mechanism by which BIG filters out A is sometimes heuristically explained on the 
basis of the behaviour of the n = 1 model in which gravity becomes higher dimensional, 
and hence weaker, at very large distances. But this model does not filter out A and the 
n > 2 model that does, has a different large distance behaviour. So let us reiterate two 
relevant aspects of the above analysis: 

1. Only a non-zero (G)(0) is relevant to filtering out A. The momentum dependence of 
(G) (that contributes to graviton decay and a varying gravitational coupling) does 
not play any role in this, in contrast to the more general "filter" mechanisms of |37| 



2. It also follows that at large distances {k 2 — > 0), the behaviour of gravity does not 
become higher dimensional, unlike the n = 1 DGP model. Hence the filtering of the 
cosmological constant cannot be attributed to such a behaviour. 

In BIG, the cosmological constant is completely filtered out at the linearized level. One 
needs to go beyond the linear analysis to make a statement about its remnant effects. In 
particular, explaining the small observed value of A (or dark energy density) is a different 
issue not addressed by BIG at this stage. 

13 Codimensions n — 1,2 give (G)(0) = oo due to an IR divergent (/-integral and the above discussion of 
filtering, formulated around flat background, does not apply. For the the n — 2 case, see |24, Eq, Gfl] . 
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4.2 On the origin of the tachyon and ghost problems 

As we have seen, a tachyonic ghost mode is encountered while analyzing BIG around flat 



background |28j. Although flat space is a valid solution to the A = model, one may still 
suspect that the tachyon and ghost problems are related to this very set up. After all, a 
generic 3-brane always has a tension which, in some sense, is the stuff that holds the brane 
together. So one always has to work with the full action, 

nbulk I abrane , abrane , abrane 

Even though a A — > limit may exist, in that limit S b £ ane = A f d A a^/—g should be 



replaced by the action for a tensionless brane [49, 50 1 rather than totally eliminated. But 
tensionless branes have special properties that make them unsuitable for our purposes and 
hence need not be considered. This line of argument tells us that the "basic" brane induced 
gravity model reviewed in section 3 does not really correspond to a brane setup, unless 
augmented by a brane tension term. Could this also provide a solution to the tachyon and 
ghost problems? This possibility is discussed below. 

One approach to addressing this question is to consider, as suggested in |3|, ^ |6], a 
large S™ ne so the matter action S^ ane can be neglected to first approximation. This 
leads to classical solutions for bulk gravity sourced by A. Around this one should study 
metric perturbations sourced by S^ ane and influenced by S"^ 716 . As pointed out above, 



the existing solutions [14, 15] are not adequate for this purpose since they do not satisfy 
the equations of motion with the correct brane source. Hence, such analysis has not yet 
been carried out (see footnote 3). 

However, the above parameter range (S^ ane >> S^ ane ) is not the only one to explore. 
In fact, we saw that the linearized analysis is already capable of exploring the model in the 
regime where S^ ane and S^ ane have comparable but small contributions. Thus, it can also 
provide some insight into the implication of the brane tension for the tachyon and ghost 
problems. 

The tachyon and ghost problems showed up in the propagator ( |3.15| ) for the s field 
which has the general structure 



3 Bk 2 - Aaj n - 2 \N(°)\Z(k) ' 

The overall sign signals a ghost and the relative sign in the denominator, a tachyon. The ap- 
pearance of A and B coefficients (associated with bulk and brane actions) shows that S(fc) 
(relevant to the tachyonic nature) arises entirely from integrating out bulk modes, while the 
rest (the —2/3 factor and hence the ghost nature) originates in the brane Einstein-Hilbert 
term. Let's emphasize that the above form is a valid way of presenting and interpreting the 
propagator as long as the B term in the action dominates over the A term. This is the phe- 
nomenologically interesting case, with a very light graviton and effectively 4-dimensional 
gravity. It is also in this range that the tachyon and ghost identifications hold . 



14 For example, in the opposite case when the A term dominates over the B term (as in string theory 
setups), the propagation is described by the S(fc) term and gravity is healthy and essentially d-dimensional. 
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The tachyon problem: Now, this is not hard to dismiss. In the previous subsection 
we saw that the "tachyonic" mass of s (more precisely, the S(0) 7^ part) was 
instrumental in "absorbing" A and keeping the brane metric flat. On the other 
hand, the classical solutions of [14, [D|] indicate that A can equally well be absorbed 



by curving directions transverse to the brane, keeping the brane metric flat. The 
comparison indicates that approximating the bulk as flat, has to force a mass on s in 
order for the brane metric to remain flat in the presence of A. In this sense, the mass 
is at least partly an artifact of the flat background approximation and is functionally 
equivalent to a "healthy" bulk curvature to be sourced by A. This demonstrates that 
the "tachyonic" mass is not an unhealthy feature. Indeed, m 2 is tachyonic only with 
respect to the ghost-like kinetic term of s, and as such has a positive contribution to 
the s field Hamiltonian ~ f(—s 2 — Vs 2 + m 2 s 2 ). 

• The ghost problem: This however is more fundamental and may not be cured even 
if the flat background is replaced by a curved one sourced by the brane tension . The 
reason is that even in such a curved background, the brane metric remains flat and 
hence its fluctuations are described by Sj^-™ 6 , expanded around flat 4-dimensional 
spacetime. This contributes a factor — (2/3)(l/Bk 2 ) to the scalar mode propagator, 
exactly as in flat background. A modified "self-energy" term, T, curve( [, arises from 
integrating out bulk modes, but now it also depends on the brane tension. As be- 
fore, small graviton masses and effectively 4-dimensional gravity require a large B 
coefficient. Thus the overall propagator again has the standard 4-dimensional tensor 
structure, corrected by self-energy terms. Now, it is well known that a propagator 
with a tensor structure associated with the 4-dimensional Einstein-Hilbert action is 
ghost free only in the massless limit (where the scalar mode ghost cancels against a 
contribution from the spin 2 mode). As soon as the mass poles are shifted away from 
zero, the scalar mode ghost remains uncanceled. 

An implication is that the ghost can be avoided only if T, curve d(0) = 0. In a flat 
space background, this would have prevented the theory from filtering out the cosmological 
constant. But in a curved background, where A is already absorbed in the background 
curvature, this could lead to a consistent theory of massless gravity with filtered out A. 
However, it is not obvious that a background sourced by the brane tension will lead to 
such a self-energy contribution. In fact the opposite seems to be the case by a continuity 
argument: In the small A limit, the expressions should reduce to what we have already 
computed, in which case S(0) 7^ 0. Hence we conclude that in the phenomenologically 
interesting parameter ranges, the scalar ghost will be a generic feature of brane induced 
gravity in any background, except those for which £ C wued(0) = 0. 

The absence of a ghost in the d = 5 DGP model is a numeric coincidence. On general 
grounds, reduction of gravity from d to 4 dimensions gives a factor d — 5 in the scalar sector. 
In the BIG setup, this results in solutions with s = 0, thus evading the ghost problem. 

4.3 Implication for couplings 



The structure of the solution (4.3) essentially corresponds to a warped metric with its 
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standard implications for the couplings. The metric in the presence of A and matter 
sources becomes, 

g ia ,(x) = (l + c)7 lia , + h$(x), (4.4) 

where c ~ A/(Au n ~ 2 ) > and off is the metric perturbation sourced by matter. Adopt- 
ing the standard convention to use t]mn as the flat spacetime metric, one defines the 
physical metric 

V (1 + c ) • ^-°) 

In terms of this, B f y/gR = (1 + c)B J ^/g'R'. For the matter action, f ^/gC(g,ip,X) = 
(1 + c) 2 f ^fg'C(g',ip',\'), where the matter fields, collectively denoted by ip, and their 
couplings A, have to be scaled appropriately for the equivalence principle to hold. This 
results in an effective B' = B/(l + c) or an effective Newton constant G' N = (1 + c)Gn, 
consistent with related observations in JL7], || || ||, £5], |2g] . Although this analysis is valid 
for small c (so that our perturbative treatment remains valid) it indicates that for a large 
A one needs a starting B larger than the observed Ml (for c > 0). This can relax the 
constraint on the bulk coupling A based on the smallness of the graviton mass by a factor 
(1 + c) which is welcome news. However, it may also introduce new hierarchy issues for 
other standard model couplings. 

4.4 Contrast with Fierz-Pauli Massive gravity 

In the ghost-free FP massive gravity theory around flat background, reviewed in a later 
section, the spin-2 graviton h 1 - has a mass m^-, whereas, to avoid the ghost, the mass tuq of 
the scalar mode s is sent to infinity. Hence for any matter source T^ u , s = 0. The caveat is 
that a constraint equation one gets for s (for tuq = oo) only implies O^s = so, in general, 
s = c\ + C2S harmonic Substituting back into the equations, gives C2 = for any source 
and c\ = — 4AG , Ar/3m 2 ! 7^ only for a cosmological constant source. To contrast the two 
theories, in BIG, s is always sourced and remains well behaved in the presence of A due 
to its own mass, but in the ghost free FP gravity, s is zero except for a A source and the 
screening parameter is the mass of spin-2 graviton m,2, which is itself not sourced by A! A 
FP theory with finite tjiq (and hence a ghost), when sourced by A, would behave similar 
to BIG for n > 2, with mo as the screening parameter. However, the mo —> 00 limit is not 
continuous. Thus, although superficially similar, BIG with n > 2 and FP massive gravity 
use different mechanisms to filter out A. 

5. Inclusion of Extrinsic Curvature Terms 

We now solve the modified equations of motion after the inclusion of the extrinsic curvature 
terms. These are ubiquitous in any brane induced gravity setup and arise at the same 



perturbative order as the induced Einstein-Hilbert term [18, IS, 52 1. The brane transverse 
fluctuations enter as new degrees of freedom. The brane tension can be included in T^ v . 
Again the analysis is performed in terms of gauge invariant variables. It turns out that the 
extrinsic curvature terms have no effect whatsoever on the solutions for the brane fields. 
While this does not cure the tachyon/ghost problem, it at least insures that 4-dimensional 
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gravity is not modified in other unwanted ways by the extrinsic curvature terms. The 
negative result is essentially due to the Z2 reflection symmetry of the background about 
the brane position. Breaking this symmetry will lead to extra non-trivial contributions 
from the extrinsic curvature terms. With this summary, the reader not interested in the 
details of the calculation can safely skip this section. 

5.1 Extrinsic curvature term in the action 



The contribution of the extrinsic curvature to the action is [18 1 

s n = c J d'a^- g (n# n A f - nt Ia «„/) , (5.1) 

where 

n% = d a d p x M - 7 ^ d x x M + r™ K d a x N d p x K . (5.2) 

Here 7 and T are the Christoffel connections corresponding to the metrics g a p{x(a)) and 
Gmn(x(ct)), respectively. For a thin brane the quadratic action for the fluctuations can be 
worked out by setting g = r\ + h, G = 77 + H and 5X M = {/**, y 1 }, 

S n = C j dx\ [~ (d a H^d a H^ - d a H ai d p H^ + 8 a H ai diH^ - d a H i/s diH als 

+ \ (diHvp&H * - diHZ&Hf) + y'di {U A H« - d a d p H a ^ +•••], (5.3) 

where the ellipsis stand for functions of h^ v — H^ u that do not contribute to the equation 
of motion for a metric variation 8h^ p = SH^ U . 

As described in section 2, for thick branes, Hmn{x((t)) are replaced by (Hmn){x\\)- 
A new feature of S(i is that it contains transverse derivatives of H restricted to the brane, 
diH\hrane = ( ~fiT~H)(x\\ > x ± = Ho)- For a thick brane, these should be replaced by (diH) = 
(■fi-rH)(x\\). The y 1 (x\\) are brane fields and do not need any "blurring". Hence, for thick 
branes, the appropriate action is, 

S n = cjdx^ (d a (H) ifi d a (H)^ - d a (H) ai d^H)^ + d a (H) ai (diH)? 

-daiH^idiH) * + ~ ((d t H) a p(d l H) a P - (d i H)%(d i H)f) 

W (u A {diH) a a - d a d p {diH) a ^ +-..]. (5.4) 

5.2 Equations of motion with extrinsic curvature contributions 

To obtain the contribution of Sn to the equations of motion for a thick brane, one considers 
variations 5y % , 6(H), 5{diH) and re-expresses the last two in terms of the bulk variation 
5H. In particular, 

8(diH) M N = / d n xj_P(xj_ - y Q ) di6H MN (x\\,x± - y Q ) 

= - d n x±d i P(x± - yo)SH M N(x\\,x± - y ) ■ (5.5) 
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Then for 5H the variation can be written in a compact from as, 



5S n = -2C d d x 
where, to linear order, 



P(x± - yo) d a N ia ? 5H ip - ^P(x ± - y ) M 1 ^ 5H a p 



(5.6) 



A^(x„) = (K P - n\ x vap) = d a d^ + - {d a (Hy p + dp(Hy a - (d*H) af) ) 



n 4 y l + d x (Hr--(d*H)$) Va p 
d a dpF* + t (d a {H% + dp{H^Y a - (d'H^^ 

d 4 f* + d x (H^) Xi - li&H^ Vafi + ^oapfts) . 

F l are gauge invariant extensions of y l (the analogues of F^ in transverse directions), 



F l = y l + (A*) + -{d l §) . 



(5.7) 



Now, from 5{S b ^ + S^g ne + S^ ane ) one obtains the thick brane metric equation of motion 
with extrinsic curvature contributions contained in the C-dependent terms, 



+P(x± - y 



±MN + d^_2( d M d N_ rj MN Dd)s 



D 



^h^ + i (d^ - fTD*) s)5^5? -Cd x Ml (S^S* + 5^5. 



iV> 



+ Cd t P(x ± - y )M t '"6?8Z = --P(x ± - y )T^S™5« • 



(5.8) 



Also, in the quadratic action Sq, the y % appear as Lagrange multipliers and their 
equation of motion is the constraint equation, 



9^m; u = ^ (n 4 # A ±A - d^H^ v + -n 4 5 



0. 



(5.9) 



But we will see that this is already contained in (|5,8| ). The expression within braces is the 
scalar curvature of the metric (77 + H)^ v . 

5.3 Solutions with extrinsic curvature contributions 

Below we show that the inclusion of C-dependent terms in ( |5.8| ) has no effect on the 
solutions for brane fields. We closely follow the steps in the C = case. 



The tjmn -trace of (5J3) gives S in terms of s and Ml v . On transforming to momentum 
space, one gets the solution 



S(k,q) 



-d 



P(M) 



A(d-l)(d-2) k 2 +q 2 \ " 2 



T^ + -Bk 2 s + i2C qi Mi 



(5.10) 
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The relevant feature of this is that, since M\ depends only on k, the C-dependent term 
is an odd function of q % and {qiM\ ) = 0. Therefore, it does not contribute to (S) which 
is, then, still given by ( |3.4| ), For later convenience, we record the analogue of (J3~5|). 



S(k,q) 



k 2 + q 2 



(S)_ 

(G) 



+ i 



2d 



A(d-l)(d-2) 



CqiN{ X 



(5.11) 



The (i,j) Components of (5.8) are not affected by the C dependent terms. Hence they 
yield again the solution ( |3.9D for H-j in terms of S. The solution ( 5.11 ) for S then gives 



the same expressions for (H- J ) and (q t q J 'H i) as in ( |3.10|) . 

Manipulations of the surface equation also lead to the same result as the C = case. 
Hence (S) and s are still related by ( 3.1 3| ) and finally the solution for s is given by ( 3.15| ). 



The (/J,,i) Components of (5.8) contain a C-dependent term and give 

d-2 k^a 1 ~ , _.C P 



H^(k,q) 



kllct S(k,q)-2i C 



Ak 2 + q 2 A 



(5.12) 



d k 2 + q 2 
However, (qiH- Lf1 ' 1 } is still given by (3.7). Hence going through the same steps as before, 



we get the same solution (|3.14| ) for F^. 

The (//, v) Components of the equation of motion in momentum space is, 

A 

~2 



(k 2 + q 2 )H^ u + d —±- (k,k u - n^ik 2 + q 2 )) S 



+P(q) 



B 



1 



k 2h U + o i k » k » ~ k ^) s ~ iCqiAf^ik) 



1, 



L [IV 



0. (5.13) 



When we solve this for H^ v and evaluate (H- 1 )^, the C-dependent term drops out (because 
of integration over an odd function of q l ). Using this to compute h^ u again gives the same 
result (|3.16 ) as the C = case. 

The y l Equation of motion (|5.9| ) follows also as a consequence of the metric equation: 
contract ( 5.13 ) with k^k v , use k^k u ii^ v = q l q J *Hh and (|3.9| ) to eliminate H^ in favour of 
S. This gives (|5.9| ) in momentum space. 

Thus, this constraint does not affect the solutions for the brane fields, but it provides 
an equation for F l as follows: Using H^ x = —H^, d^d^H- 
rewrite these in terms of S, one gets 



' \IV 



&8>H± and Q to 



k^N;„ = --(d-l)(d-5)k 2 {q*S) 



0. 



(5.14) 



Ignoring the 4-dimensional harmonic solutions, this gives (q l S) = 0, or (d l S) = 0. Using 
the solution for S in ( |5.10| ), this becomes (with no sum over j), 



(Q J S) 



which holds only if, 



2dC 



A(d-l)(d-2) 



d n q P 2 (q)W? 
k 2 + q 2 



*l 



0, 



M{\k) = 3k 2 Fi + 3i{ qi H ±ji ) - ihq j H^) + A(q j S) = , 



(5.15) 



(5.16) 
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Again, expressing H^j in terms of 5 gives 



^P w+5 ^Oh°- 



But ( 5.11| ) implies (qi S/(k 2 + q 2 )) oc A/j[ = 0, so we finally get the F % equations, 

□ 4 F J ' = 0. (5.17) 

Hence, up to harmonic functions, F l = 0. This is a consequence of the reflection symmetry 
of the set up in directions transverse to the brane. A background that breaks this symmetry 
can lead to non-trivial F l . Note that for the DGP model with d = 5, the constraint equation 
(5.14) is trivial and does not imply an equation for F 5 . 

Thus the presence of the extrinsic curvature related corrections has no effect on BIG at 
the linearized level, in particular not on the ghost problem. While this is discouraging, the 
positive aspect is that such terms (that on general grounds are ubiquitous in all braneworld 
models) do not cause large deviations from standard gravity. 

6. Revisiting Fierz-Pauli Massive Gravity 

In the remaining part of this paper, we turn to a comparison of brane induced gravity for 
n > 2, with massive versions of Einstein-Hilbert gravity in 4 dimensions. In this section we 
start with a review Fierz-Pauli massive gravity theory [p9| |30| (for reviews, see [^, [5j]). 
The theory is then rewritten in a way that facilitates comparison with the 4-dimensional 
effective action for BIG, to be derived in the next section. In fact, the construction presented 
here for the FP theory enables us to make sense of the BIG effective action. 

6.1 Fierz-Pauli massive gravity in gauge invariant variables 

In standard 4-dimensional Einstein-Hilbert gravity, metric fluctuations h^ u in a flat back- 
ground can be made massive by introducing mass parameters, say, a and b. The com- 
bination a + b = results in the ghost free Fierz-Pauli massive gravity theory |2S, |30| . 
Since we are interested in the origin of ghost in brane induced gravity, we keep a and b 
unconstrained and use the term "Fierz-Pauli theory" in this generalized sense. 

The gauge invariance of linearized gravity, broken by the mass terms, is restored by the 
Stiickelberg technique. This amounts to performing a gauge transformation and retaining 
the gauge parameters as new fields. The resulting action is, 

Sfp = S EH [h] ~jjd 4 x [o( V + 23 (M /,)) 2 + b(h» + 2djn 2 ] • (6-1) 

Seh [h] is read off from ( |2.6| ) for d = 4. The / M are Stiickelberg fields transforming as <5/ M = 
A^. Combined with bh^ v = —2duX l/ \, it keeps Sfp invariant. The /^ can be interpreted 



as the Goldstone fields corresponding to the broken symmetry [51|. This conveys a sense 
of graviton masses arising due to spontaneous symmetry breaking, even in the absence of 
a more detailed fundamental mechanism. 
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The theory can be rewritten in terms of h and s of 
variables, 



and the new gauge invariant 



(6.2) 



Note that these are the same as the F^ in ( 2.17 ) once we identify the Stiickelberg field as 
f/i = ffi J r (Aft) + id M (($}), hence the similar notation. The Fierz-Pauli action ( |6.1| ) in 
terms of the gauge invariant variables becomes, 



S FP = -| j (fx\h^{-U + a)h^ u + s{^U + ^ + b)s + A{^ + b)d ■ F s 
+ a {d„F v - 8 U F^) 2 + 4 (a + b) (8 ■ F) 2 ) 
The equation of motion for F^ , 

a DF M + (a + 2b)d ll {d ■ F) + (~ + 6)d„s = , 



has the purely longitudinal solution, 



2 -4-fo 



2(a + 6) 



(6.3) 



(6.4) 



(6.5) 



Then, the L„ equation (obtained after expressing variations 5k 1 - and <5s in terms of 5h 
using the projection operators in Appendix A), 



(-□ + a)fc AM , + -(»7,H, 



a - 



has the massive graviton solutions, 

1 1 ~ 



h 



(41/ 



(|n + | + 6)* + 2(| + 6)a-F) -■/;,„/£. u>.o) 



(6.7) 



corresponding to masses wi?± = a and m^ = — ^ ^3 (a + 46). The full metric is obviously, 



/;. 



/"' 



^ + 7 V^ + d V F " + ^ F M . 



with the gauge choice /„ = 0. In the following paragraphs, we comment on the ghost free 
theory before returning to the general case in the next subsection. 

The massive FP gravity theory is, in general, not ghost free: A standard analysis of 
the gauge invariant amplitude ( |3,17 ) for interaction between sources T^ v and T' shows 



that for a /0, the s field is always a ghost , as shown, for example, in [29, 31, 44 1 or 



D The presence of the ghost for a + b 7^ can already be understood at the level of action (3.3) in 



terms of the Stiickelberg fields: The vector field F^ = F^ + d^ipF has an unhealthy kinetic term because 
(a + fe)(<9 M F M ) 2 leads to a 4-derivative term for the longitudinal mode 4>f, resulting in ghost non-decoupling. 
For a + b — 0, the 4-derivative term is eliminated. This provides an intuitive way of detecting the massive 
gravity ghost. 
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more explicitly in [45]. The only way to get rid of the ghost is to set a + b = which gives 
m s = oo and s = 0. The outcome is the ghost-free Fierz-Pauli massive gravity with 16 , 

1 i 1 m± . ^ k^k u 



An alternative solution to the ghost problem is to regard gravity as an effective theory and 
give the ghost a mass above the scale of validity of the effective theory by setting a + b = 5 
where 5 is positive and very small. Then the inaccessible ghost pole becomes physically 
irrelevant | |55| |. 

Let us discuss an issue raised in section 4.4. While the ghost free theory has s = for 
a generic source, there is a caveat when T^„ corresponds to a cosmological constant: When 
a + b = 0, ( |6.4[) implies D^s = and its solution s = c\ -\- C2S harmomc replaces the expression 
for s in terms of d ■ F obtained from (|6.5| ). In general c\_i = by boundary conditions, but 
for a cosmological constant source, the h^ u equation gives s = c\ = — 4A/3aB. In section 
4.4, this behaviour was contrasted to that of brane induced gravity. For generic a and b 
the two theories respond to the cosmological constant in the same way. 

6.2 A 1-parameter family of actions for FP gravity 

From the outset it was clear that massive gravitons and ghosts were common features in 
both Fierz-Pauli theory (for a + b ^ 0) and BIG (for n > 2) at the linearized level. Since 
FP gravity is simpler, it is natural to use it as a prototype to model some aspects of BIG, 
as in J38|, p^| . But how much can the similarities be stretched? To quantify this, below we 
rewrite the Fierz-Pauli action in a way that helps us compare with, and make sense of, the 
BIG 4-dimensional effective action, to be derived in the next section. 

As far as classical solutions are concerned, the action ( |6.3| ) is not unique. To find other 
equivalent actions, use the relation between s and d ■ F obtained from ([0]) to convert a 
part of the s 2 term into (d ■ F) 2 and sd • F terms. This gives, 

S FP = "f / d 4 x \h^ v {-U + a)h^ u + ^sDs + (| + b + w^j ss 

+ 4(| + b + w 2 )d -Fs + a (d^F v - d v F^) 2 + 4 (a + b + w 3 ) (d^) 2 ) ] , (6.8) 

where, wis 2 + 4ui2d ■ Fs + 4^3(9 • F) 2 = should hold on the classical solution. The tD's 
are further constrained by requiring that the above action leads to the correct equations 
of motion for F^ and h^ iu . The outcome is, 

a+b _ _ a+b _ 

a/4 + b a/4 + b 



16 The origin of vDVZ J3l| , [32] , p3| discontinuity, that arises in the massless limit of ghost-free FP theory, 
is manifest in this formulation: Sending m h ± — > 0, keeping m s — oo fixed, does not recover the standard 



massless solutions (3.18) as s remains zero and the theory does not couple to Tff. The difference between 
the two massless theories does not affect light for which T^ = 0. But matter is affected differently and the 
definitions of the Newton constant in the two theories differ by a factor of 3/4. This shows up in comparing 
the bending of light in the two cases. The discontinuity is avoided in the non-linear theory through the 
Vainshtein effect [ElEct 
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Thus we get a 1-parameter family of actions, in terms of gauge invariant variables, pro- 
ducing the same equation of motion as the Fierz-Pauli theory. 

A point to note is that, except for w\ = 0, these equivalent actions will look non-local 
when expressed in terms of hn V , with D _1 -type non-localities arising from the projection 
operators (A..1) in Appendix A. 

The particular choice W3 = —(a + b) (implying wi = —(a/4 + 6)) decouples gravity 
from F^ which become a free Abelian gauge field. The resulting action involves only h 
and s (corresponding to integrating out F^), 

S F P = -^j <?x [h^i-U + «)/,£, + !*(□+ ^f ) «] • (6-10) 



Note that, although the solutions ( p.7j ) can be readily obtained from here, this choice of 
W3 is a singular limit of ( |6,8| ) in the sense that d ■ F becomes undetermined 17 . 

When re-expressed in terms of huv, using projection operators of Appendix A, this 
action remains manifestly gauge invariant, with no dependence on Stiickelberg fields and 
no sign of symmetry breaking. But, as a trade off, it will contain D -1 -type non-localities 
arising from projection operators. In this example, we know that these non-localities arise 
as a result of integrating out the F^ and indicate the existence of a local formulation (( |6.1| 
or (B^)) in which Stiickelberg fields are needed to restore gauge invariance. The two cases 
are connected by the interpolating action ( |6.8[ ), This construction is useful for obtaining 
the 4-d effective action in brane induced gravity where the above steps are traversed in the 
reverse order. 

7. 4-d Effective Action in BIG and Broken Gauge Invariance 

In order to better understand the 4-dimensional structure of BIG, we need an effective 
4-dimensional description of it. In this section we work out such an effective action by 
integrating out parts of the fields with support away form the brane as well as the extra 
components of the metric. Then we construct an action containing the Stiickelberg-like 
fields Fp, based on our understanding of massive gravity in the previous section. 

This has implications for the origin of ghost and for the realization of the brane general 



covariance as a spontaneously broken gauge symmetry: The F^ or /„ ( 2.17 ) of BIG are 
closely related to Stiickelberg fields / M in Fierz-Pauli theory (more precisely, / M ~ f^ + 
(Ajj) + ^<9 M ((<I>))). Now, fn restore gauge invariances broken by the FP mass terms and, 
hence, can be interpreted as the Goldstone fields associated with the broken symmetry 
1 51]. The appearance of their analogues /^ in BIG is thus a manifestation of spontaneous 



breakdown of 4-dimensional general covariance by the bulk-brane setup. In this sense, BIG 
provides a realization of the gravitational Higgs mechanism. 

7.1 4-dimensional effective action in Brane Induced Gravity 

In this subsection we derive a 4-dimensional effective action for BIG by integrating out all 
bulk related modes in the A and B terms of (|2.1| ). The bulk field Hmn(x\\,x± — yo) has 



This is obvious from the _F M equation of motion for (|6.8|) which gives, 2(a+b+W3)d-F = — (a/4+b+W2)s 
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a part (H- l )mn{x\\), given by ( 2.10J ), that has support only on the brane. The remaining 
part AHmn(x\\,x± — yo) has no support on the brane, i.e., (AHmn) = 0, and is given by 
the decomposition, 

Hmn = AHmn H ^ (H)mn, (7.1) 



Z 



where, Z = J d n x± P 2 = (2ir) n J d n qP 2 . In terms of momentum space variables, 



H M N(k,q) = AH M N(k,q) + 



Z 



(H) MN (k) 



(7.2) 



To determine the effective action on the brane, we eliminate, first, AHmn by using its 
equation of motion and then solve for {H)ij and (H)^j using the corresponding equations. 
This leaves us with an effective action for (H)^ u or, equivalently, an action for the 4- 
dimensional brane variables h^^s, and F^. 

To integrate out AHmn, start with the bulk action ( |2.6D , with Hmn given by the above 
decomposition, and a Lagrange multiplier term implementing the condition (AH)mn = 0, 



S% k + J d d x(PAH MN ) L MN ( 



(7-3) 



To insure gauge invariance, L MN must satisfy 8m{PL mn ) = 0. The momentum space 
equation of motion for the variation 5Hmn = 5(AHmn) (with no brane support) is, 



A 
~2 



2£-±MN | « ~~ 2 w v 
P tl -\ 



(p M p» - 7/ M V) S 



L Mn (k)P(q) = 0. 



(7.4) 



This can be easily solved for S and H in terms of L = L¥, and L . Blurring 



the solutions using ( 2.11| ) then determines L and L in terms of (S) and (H ) . 
Substituting back in the solutions for S and H ±AIN gives (|3.5|) and , 



# 



±AfiV 



P 



■(^) 



A/ AT 



+ 



d 



\MN 



d 



k 2 + g 2 (G) 
where we have introduced (as a generalization of ( p.8[ )) 

p M p Np2 {q) 



+ 



p M p N 



(G) k 2 + q 2 



(S) 



(7.5) 



(?) 



MN 



(2ir) n / rf 



(fc 2 +<7 



2\2 



(7.6) 



Substituting this back in S^ gives a 4-dimensional action for (H ) and (S) (the 
contribution of the Lagrange multiplier term vanishes) , 



q 



bulk,eff 
EH 



CM 



,a 



d 4 k- 



i 



<^> mjv <^>a™ - ^'^''h sHsy 



d 2 

d-2 (g)MN r r T ± ,m x 



(G) 



(H^) MJ}I {S)* . (7.7) 



From this we can further eliminate (H }%j and (H )%n by expressing them in terms of (S) 
using the solutions of the (i,j) and (/U,j) components of (3.2) along with ( |3.5| ). Then, we 
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express (H )fj, v an d (S) in terms of hh v , s and F^ using the surface equation (|2.16 ) and its 



consequences. Note that the equations used so far also determine F^ in terms of J ( 3.14| ), 
but leave s and h 1 - undetermined. Finally, adding 5'^g ne [/i] leads to the 4 -dimensional 
effective action, 



S eff [h} = {2^j l Vk 



#+*®-^2+\{*+^®' 



(7.8) 
This gives the correct equation of motion with the correct solutions ( 3.15|) ,( p.l6| ). However, 
it is not yet the final form of the effective action, as explained below. 



The above action is the analogue of the Fierz-Pauli action in the form ( 6.10 ) which, as 
argued there, is not its most natural {i.e., most local) form. In fact, the /i^-terms in the 
two actions can be mapped by replacing the FP mass parameter a as, 

a^a(k) = -^(G)(k). (7.9) 

Naively, the s-terms too can be mapped by replacing the FP b parameter by a(k)/(d — 6). 
But that is not the correct comparison since it does not map the F^ solution of FP theory 
(|6.5|) to the F^ solution of BIG (3.14) (simply because the quantity U does not appear in 



the map). In spite of this, the similarity between the BIG action (|7.8| ) and the FP action 
in the form ( |6.1C ), teaches us the following: 



Recall that in the FP case, ( |6.10| ) was a singular limit of a 1-parameter family of actions 
( |6.8j ), a generic member of which contained the FP Stiickelberg fields Fa. The w\ = 
member of this family was the original FP action (|6.1| ) which was local when written in 
terms of h^ and f^. However, the specific w\ choice that led to ( ft.lOQ decoupled F^, 
so the corresponding solution could not be reproduced. Also this action had l/A; 2 -type 
non- localities when expressed in terms of hnv- 

To compare, this also is the structure of the BIG effective action ( |7.8| ). It too does 
not contain the worldvolume fields Fa and, when expressed in terms of hnvi W1 ^ contains 



l/fc 2 -type non-localities arising from the projection operators ( A.l ) 18 . We take this as 
an indication that ( |7.8| ), too, is the singular limit of a 1-parameter family of actions that 
generically contain the i 7 ^, beside the metric. 

The fact that the replacements for a and b described above cannot map the F^ solution 
to the Ffj, solution means that the 1-parameter family to which the BIG effective action 
Q7.8D belongs cannot be obtained from the corresponding construct ( |6.8| ) in FP theory. In 
this sense, the two 1-parameter families are not equivalent. The explicit construction of 
the BIG 1-parameter family is performed below and is then compared to the FP theory 
for a different choice of b. 

7.2 A 1-parameter family of 4-d BIG effective actions 

To reinstate the F^ into the BIG effective action, let's get back to equation Q7.7|). We again 
eliminate the bulk tensors in favour of h^, s and F^, but now retain F^. Since the relation 



IX, 



Of course, the action (7.8) also contains non-localities through the fc-dependence of k (G). But these 



arise as a result of integrating out bulk modes and are not related to Stiickelberg fields. 
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between s and F^ is already implicit in the equations employed, the outcome is not unique 
and depends on the steps followed. Depending on the details, the resulting action may not 
even lead to the correct equations of motion for h^ u and F^. But this can be corrected. 
Indeed, a candidate action that one directly obtains in this way is 19 , 



with, 



o 



bulk.eff 
EH 



(M 



A 



d*k 



L(G) M (G) 



+ JL ( k 2 F 2 + ( fe . f) 2 ) + i -S- ( k ■ F*I 



(G) 



(G) 



k-FT 



(7.10) 



U-i 



(d-2) 



16 (d-l) 2 (d-5) 2 



2(d-4)(2d-5)-9- 



d-1 



3 (d - 7) 



(G) 



U- 



2 — T 



(d-2) 



4 (d-l)(d-5) 



2(d-4) -3-3-^- 

(G) 



(7.11) 



This does not reproduce the correct equations of motion (after adding <S^g ne [/i]), indicating 
that the ambiguity between s and d ■ F terms has not yet been properly resolved. To cure 
this, one can use the relation between s and d ■ F in ( 333 ) to write an appropriate part of 
the s*s term in terms of k ■ F. As a shortcut, a simple inspection shows that the correct 
solution F^ = —iU(k fM /Ak 2 )'s can be obtained if TJi is, somehow, replaced by U given by 



(3.14). Hence, the trick is to set \J 2 = U + 5U2 and then convert the coefficient of the 
8U2 piece to ss* form using ( |3.14j ). This leads to an effective action containing F^ that 
reproduces the correct equations of motion, 



S eff [h,F]=S%r[h] + (2n)^ 



d"k 



r _L h ±n» h ±* + _L 

(G) ' 4(G) 



3 d 



+ 



(G) 



k 2 F 2 + 



FY +i 



u 

W) 



4 d-5 
k-F*s-k-Fs ' 



+ U' 



ss 



(7.12) 



But, as experience with FP theory shows, this action may not be unique in this respect. 
Consider converting a part of the ss* term into (dF) 2 and dFs terms, thereby, shifting the 
corresponding coefficients by u>x, 0J2 and W3 to (for convenience we use the notation Q7.9|)), 



S eff [h,F] = -(2vr) 4 |/d 4 A; \(k 2 + a)h^ v h 



"fJ,U 



w 



5(A^ + ¥+«i) " 



i{U + w 2 )a (k-F*s- k ■Fs*)+ a{k^F v - KFrf + a{A + w 3 ) (k ■ F) 2 . (7.13) 



0. 



Clearly, the modification is constrained by w\ss* + iw2(k- F*s — k ■ Fs*) +Ws(k- F) 2 
Then the equations of motion reproduce the correct solutions provided, 



w$ 



U 



w 2 , 



U>2 



u 



W\ 



(7.14) 



19 We have rewritten k^k" (H ± }^ lv {S) using the surface equation (2.16). Instead, writing k^k" {H ) IJbV — 
(qiqjH ±l:i } and using the second equation in (|3.1C|) would have given a different result. 
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Hence, we have a 1-parameter family of correct BIG effective actions. 
The F^ equation of motion, 

2k 2 F u + (2 + w 3 ) k a (k ■ F) + i(U + w 2 )k^s = , 



gives the correct F^ solution ( p. 14 ), except for wi = — U (implying ^3 = — 4) which leads 



again to the action (|7.8[). For this special value, d^ l (d ll F I/ — d l/ F ll ) = and F^ is pure gauge, 
but otherwise undetermined. Hence, as expected, (|7.8f) arises as the singular limit of the 



family ( 7.13 ). The the complete effective action, capable of determining all worldvolume 
fields h^, s and F, should be a non-singular element of this family. 

The above discussion also highlights a structural similarity between the BIG and FP 



actions in the forms ( 7.13 ) and (6.8) which are inequivalent otherwise: The specific choices 
of the parameters Wi and w-i that decouple F u and F u , respectively, from gravity, also 
convert the F J s into free Maxwell fields. This, in particular, decouples the longitudinal 
mode of the vector field which is ghost-like, thereby, transferring the ghost entirely to 
the gravity sector. In this respect, the ghosts in BIG (for n > 1) and in FP theory (for 
a + b 7^ 0) have similar origins, both being related to the survival of the gauge dependent 
components of the brane metric (contained in F^). 

7.3 Spontaneous breaking of 4-d gauge invariance 

The discussion in the previous section shows that the correct BIG effective action is a non- 



singular member of the family (7.13). The presence of F^ in this action has the outcome that 
gauge dependent components of the metric, a^ and <fi ( |2.8| ), do not drop out of the action and 
gauge invariance is maintained only through the compensating transformations of fa- Thus, 
the 4-dimensional gauge invariance of the theory, i.e., symmetry under 5h uu = —2dt^X v \, 
which is manifest in the starting BIG action, is broken spontaneously in the effective 4- 
dimensional theory, as signaled by the need for compensating transformations. 

This effect can be traced to the fact that the process of integrating out non-brane 
components of graviton, by partially solving the equations of motion, requires eliminating 
gauge degrees of freedom through either explicit gauge fixing, or the use of gauge invariant 



variables. In our approach, this is technically related to the necessity of using ( [2.16 ) in 



relating gauge invariant bulk and brane variables. In other approaches, the same issue will 
arise though it may be obscured by an explicit gauge fixing. 

To make the spontaneously broken nature of the symmetry very explicit, we can write 
the effective action in terms of the gauge non- invariant variables h^ y and /„. A familiar 
form is obtained by matching some terms of the BIG effective action ( |7.13 ) (by choosing 



an appropriate w\) to the starting form of the FP action ( |6.3| ). For example, we can map 
the coefficients of the s 2 and (dF)s terms to FP theory by choosing 

3 d-2 U U 2 

Wl ~ ~16d-5 ^U + 2> ~ T ' 

Then, in terms of 

a(k) = ~{G)(k) , b{k) °< :if] -' 2 ' 



B x /v " v ' 4 V2 d-5 U + 2 
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the BIG effective action takes a form very similar to the FP action, differing from it only 
in the (dF) 2 term. Finally, in terms huv and /^, it becomes, 



S eff [h,f] = S EH [h}-^Jd 4 k 



a(k) V + 2* k^L) + b{k) h» + % W 



2U + 1 



U 



1~ 2ki"k 



V 



2b(k) k-f + -h -5-/1 



6 3 k 2 



ixu 



(7.15) 



The first line is an exact map from the FP theory. The second line shows the structural 
difference between massive FP gravity and the BIG effective action for n > 1. 

Finally, one may comment on the similarity of brane induced gravity with the gravita- 
tional Higgs mechanism of (|35|, J36|). In these works, in order to break general covariance, 
one introduces 4 scalars fields </>^(cr) = cr^ + d^ , with a Minkowski metric on the field space 
(er^ denote 4 dimensional spacetime coordinates). These fields are very similar to our brane 
embedding coordinates, x IM (a) = a^ + f 1 * which naturally come with a Minkowski metric 
(the transverse fluctuations y l are not relevant). In this sense, BIG provides a natural 
realization of Higgs mechanism in gravity. 

Note that the presence of extra dimensions is crucial for this similarity to hold. The 
reason is that, as discussed earlier, the /** are non-trivial only for n > 0. It is interesting 
to explore the connection between these two theories further. 

8. Conclusions 

Our results are summarized and discussed in detail in section 1. Here we summarize 
the main conclusions. The existing classical solutions for bulk gravity sourced by a 3- 
brane tension are not adequate for studying metric fluctuations on the brane. Hence we 
use the flat background approximation but argue that the results also shed light on the 
curved background analysis. In the flat approximation, brane induced gravity filters out 
a cosmological constant using the mass of the unhealthy scalar mode. This mass is partly 
an artifact of the flat approximation and there is no real tachyon instability (in the sense 
of an unbounded Hamiltonian) associated with it. 

The presence of ghost in BIG is related to the phenomenologically interesting values 
of parameters in the theory, corresponding to a very low graviton mass and effectively 
4-dimensional gravity on the brane (for example, there is no problem with the natural 
parameter values in the string theory D-brane setups) . For these parameter values avoiding 
the ghost is difficult even in a curved background sourced by A since our analysis can already 
probe such backgrounds for small values of A and the problem is very similar to ghosts in 
ordinary massive gravity. If this turns out to be the case, then BIG cannot be regarded as 
a viable solution to the cosmological constant problem. 

The only healthy exception would be a background that absorbs the cosmological 
constant by curving directions transverse to the brane while leading to massless gravity 
on the brane (T, curvec i(0) =0), if that is possible at all. This would avoid the ghost issue 
and lead to a theory with filtered out A and modified massless gravity on the brane. The 
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modifications to gravity can be computed using the methods of section 3.3 and wili have 
observable consequences. However, verifying this possibility requires constructing brane 
classical solutions that satisfy the equations of motion with the correct source terms. As 
indicated earlier, finding such solutions is a generic problem in brane physics. At the 
moment, we cannot say more on this possibility (however, see see footnote 3). 

Another conclusion is that, using gauge invariant variables, it is possible to keep track 
of symmetries and explicitly show that, in the effective 4-dimensional gravity theory, gauge 
dependent components of the metric do not decouple. This indicates a spontaneous break- 
ing of 4-dimensional general covariance as a consequence of the bulk-brane setup. This 
setup generates the extra gauge invariant modes that contain the Stiickelberg fields as- 
sociated with the spontaneous breakdown of symmetry. It provides a realization of the 
gravitational Higgs mechanism. 
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A. Projection operators for gauge invariant variables 

For a symmetric tensor Hmn, its transverse part H' MN satisfying d M H' MN = is given by 

H MN = H MN - Dj\d M d L H LN + d N d L H ML - a d 1 d M d N d L d K H LK ) 

Then, the components in the decomposition ( |2.5| ) are given by (with H = H^), 

± , 1 ( d M d N , d d L d K H LK 
H M n -Hmn-j^ (vmn - ~^-)H , S - — (H ) , 



a X txLjj o d L 8 K H LK . d ( d L d K H LK ^ 1 H 



(A.l) 



The gauge transformation 5Hmn = — 2d( M £/v) results from the variation of Am + \3m&- 
H MN and S are invariant except under a restricted class of transformation with Od£ = 0. 
These become relevant when D^-Hm^ = which is not the case in brane induced gravity. 
The energy-momentum tensor Tmn can be decomposed in a similar way, 

Tmn = T MN + BmTn + OnTm + OmQnT^, + -tjmnTs 

The components in this expansion are given by expressions analogous to those for Hmn in 
terms of Tmn and a T' MN . 

Generically, OjTmn ^ 0. Then 8 m Tmn = implies T' MN = Tmn and other sim- 
plifications. However, if O^Tmn = 0, then D^" d L d Tm / and one should not use 
the simplified equations. An example is the cosmological constant source T MN = tjmn^- 
Then, the simplified equations do not hold, while using the complete projection operators 
one gets, T^~ N = T M = T, A = and T A = dk. 
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